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Abstract 

We extend the standard model to a scalar-assisted vector-like fermion model to realize electroweak 
baryogenesis. The extended Cabbibo-Kobayashi-Maskawa matrix, due to the mixing among the vector¬ 
like quark and the standard model qnarks, provides additional sources of the CP violation. Together 
with the enhancement from large vector-like quark mass, a large enough baryon-to-photon ratio could 
be obtained. The strongly hrst-order phase transition could be realized via the potential barrier which 
separate the broken minimum and the symmetric minimum in the scalar potential. We investigate in 
detail the one-loop temperature-dependent effective potential, and perform a random parameter scan to 
study the allowed parameter region that satisfies the strongly first order phase transition criteria Vc > D- 
Several distinct patterns of phase transition are classified and discussed. Among these patterns, large 
trilinear mass term between the Higgs boson and the scalar is preferred, for it controls the width of the 
potential barrier. Our results indicate large quartic scalar couplings, and moderate mixing angle between 
the Higgs boson and the new scalar. This parameter region could be further explored at the Run 2 LHC. 
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1 Introduction 


The baryonic matter that remains after the baryon-antibaryon annihilation, makes up around 5% of the 
total energy density of the universe. It is puzzling that the universe does not have equal amounts of matter 
and antimatter. We can characterize the asymmetry between matter and antimatter in terms of the baryon- 
to-photon ratio 


riB 

V= —, 


( 1 . 1 ) 


rij 

where ns = nj, — is the difference between the number density of baryons and antibaryons, and is 
the number density of photon. The rij is introduced to prevent the parameters r] from diluting during 
the expansion of the universe after nucleosynthesis. The baryonic matter desity ub sd, present time has 
been consistently measured by the big bang nucleosynthesis and the fluctuations of the cosmic microwave 
background. The Planck result for the cosmological density parameter [1] 


flBh^ = 0.02226 ± 0.00016, 


( 1 . 2 ) 


translates to the baryon-to-photon ratio 

77 = (6.05 ±0.07) X 10"^°. (1.3) 


Explaining the observed baryon asymmetry has been one of the greatest challenges of particle physics 
and cosmology. As the entropy production during inflation could greatly dilute and thus wash out any 
existing baryon asymmetry, it is reasonable to assume a zero baryon number density after the inflation. 
Later, the asymmetry is generated dynamically through the so-called ” baryogenesis”. It has been suggested 
by Sakharov [2] long time ago that the general baryogenesis has three necessary conditions: baryon number 
violation, sufficient C and CP violation, and departure from thermal equilibrium. Hence, we look forward to 
a mechanism in which these three conditions are satisfied and could provide the observed baryon asymmetry. 

Electroweak baryogenesis [3-5] (EWBG) offers a theoretically attractive and experimentally testable 
mechanism to realize baryogenesis. The great attraction of this mechanism is that the baryogenesis took 
place at or near the electroweak scale, suggesting that it might be probed in the near future by experiments 
at the accelerators. The EWBG proceeds as follows (see [6-11] for reviews). At temperatures far above 
the electroweak scale, the electroweak symmetry is manifest, which implied a high sphaleron rate that 
preserves baryon symmetry in thermal equilibrium. As the universe cools down to near the electroweak 
phase transition scale, bubbles of the symmetry-broken vacuum began to emerge and grow. CP violating 
processes involving the electroweak sector were triggered at the expanding wall of the bubbles, leaving baryons 
inside the bubbles and antibaryons outside. Through the rapid sphaleron transitions in the unbroken phase, 
the excess of antibaryons are washed out. Meanwhile, if the sphaleron rate in the broken phase could be 
suppressed enough, the excess of baryons inside the bubbles could survive. We can easily realize a Boltzmann 
suppression of the sphaleron rate, because the sphaleron has an excitation energy Esph that is related to the 
Higgs vacuum expectation value (vev) v. It has been shown [12] that the suppression is strong enough when 

Esph/Tc > 45, (1.4) 

which serves as the condition for a strong first-order phase transition (SFOPT) in the context of electroweak 
baryogenesis. 
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The standard model (SM) contains all the necessary ingredients to realize electroweak baryogenesis: 
baryon number is violated by sphaleron processes; CP violation comes from the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix; departure from equilibrium is realized by the bubble nucleation and expansion during the 
first-order electroweak phase transition (EWPT). However, given the observed Higgs boson mass Mh = 125 
GeV, the EWPT is not strong enough to suppress the sphaleron rate inside the bubbles [13]. Also, the CP 
violation in the CKM matrix is not large enough to generate the expected asymmetry. Therefore a successful 
electroweak baryogenesis needs new physics beyond the Standard Model. The new physics should provide 
new sources of CP violation that can be manifested by the advancing bubble walls, and also provide strong 
enough first order EWPT. Both conditions require the existence of new physics at around the electroweak 
scale that directly couples to the SM Higgs sector. A simple and economic way to realize the strong first order 
EWPT is to add a new scalar which couples to the Higgs boson, such as the singlet extended standard model, 
etc [14-22]. Moreover, if the scalar is a real singlet [18,22], the cubic terms could exist in the potential at 
tree-level, and therefore the phase transition gets stronger without the need of the thermally induced barrier. 

We consider the electroweak baryogenesis in a scalar-assisted vectorlike fermion model [23], in which a 
singlet scalar and vectorlike fermions are added to the SM particle content. Originally, the model is motivated 
by the possible instability of the vacuum structure in the vector-like fermion model. The singlet scalar is 
added to the scalar sector and couples to the vector-like fermion. This model solves the vacuum stability 
problem in vectorlike fermion model and possible perturbativity issues in singlet scalar extended standard 
model. Recently this model attracts lots of attentions because it could naturally explain the diphoton excess 
observed at both ATLAS and CMS [24, 25] . The diphoton signature of this model and its extensions has 
been considered in Ref. [26]. Due to constraints from other WW and ZZ channels, a 750 GeV scalar singlet 
could accommodate the observed diphoton excess more readily than the SU(2)l x U{1)y scalar multiplets. 
The diphoton signature is produced via the gluon fusion and subsequent diphoton decay with vector-like 
fermion running in the loop. 

In this work, we consider that this model realize the ectroweak baryogenesis. The vector-like fermion 
mixes with the SM quarks, extending the 3x3 CKM matrix to a 3 x 4 matrix, which provides additional 
sources of CP violation. Due to the coupling between Higgs and the new scalar, the phase transition happens 
in an extended scalar space, which leads to more possibilities on phase transition. We will discuss the scalar 
potential in detail, perform numerial calculations, and investigate how the extended scalar sector provides us 
the SFOPT. Furthermore, we classify the phase transition patterns and explore the parameter preferences 
in each pattern using the shape of the derivatives of the scalar potential. Finally, we explore the discovery 
potential of the parameter space favored by SFOPT at the LHC. 

The organization of presentation is as follows. We begin with the description of the new physics model. 
In Sec. 3, we discuss the CP violation in this model. In Sec. 4, we present the effective potential and the 
shape of the scalar potential. In Sec. 5, we discuss the phase transition pattern and explain our numerical 
results. In sec. 6, the LHC discovery potential of the favored parameter region is discussed. We then make 
our conclusion. 
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2 Scalar-assisted Vectorlike Fermion Model 


In our setup, we consider an extension of the SM in which a vector-like fermion U and a real singlet scalar s 
are added to the particle content [23]. The vector-like fermion U transforms as (3,1)2 under the SM gauge 
symmetry SU{3)c x SU(2)l x U{1)y- Due to the same quantum number, its right-handed component mixes 
with the SM right-handed up-type quarks. As is known [23], the vector-like fermion model encounters the 
vacuum instability problem. To have a stabilized scalar potential, a singlet scalar s is introduced. The scalar 
mixes with the SM Higgs boson and couples to the vector-like fermion. Here, we assume no Z 2 symmetry for 
the new scalar, so that it has a non-zero vev in general. We refer this model as the scalar-assisted vector-like 
fermion model. Let us discuss the quark sector and the scalar sector in detail. 

In the quark sector, a SM family contains a doublet Ql = {uL,dL)^ ^ (3,2)i/6 and two singlets un ~ 
(3,1)273 and dfi ~ (3,l)_i/3 that couple to each other via the Higgs doublet $ ^ (1, 2)1/2 • Due to the 
flavor experiment constraints, we only allow the new vector-like fermion to have significant mixing with the 
third generation quarks, so it’s also called a top partner. Vector-like top partner is well motivated in the 
little Higgs [27], composite Higgs [28], and extra dimension models [29], etc. Let us first write down the 
Lagrangian for one flavor mixing between the new fermion and the third generation quarks, and then extend 
to the three flavor mixing. We could write down the following new Yukawa couplings: 

^Yuk = ytQL 3 ^ — y'Qm^ — ygsU lUr — MUrUr + h.c., (2.1) 

where Ql 3 and urs are the left-handed quark doublet and the right-handed up-type quark in the third 
generation. The vacuum expectation values of the two scalars are denoted as 


V={4>), U={s). 


The mass term becomes 


r — — 

■^mass — 


(n„ C.) (7 


vy 


Ur 


+ h.c.. 


0 M -I- ysu) \Ur^ 

To get the mass eigenstates (t,T), we diagonalize the fermion mass matrix 

Ivyt vy’ \ ( cosBl sind^N /to* 0 \ /cos0/j -sinOR^ 

\ 0 M + yguJ \ —sin^i cos9lJ V 0 mxJ \sin0fl cos9r , 


Note that the two mixing angles are not independent parameters. 


( 2 . 2 ) 


(2.3) 


(2.4) 


tandjj =- tanOx. (2.5) 

rriT 

Despite the tight constraints on the flavor mixing between the new vector-like fermion and the first two 
generations, these mixings are still essential for the new CP violation. If we consider the three families of 
the quarks in the SM, the Yukawa couplings yt and y' in Eq. 2.1 becomes matrix and vector Y/ in the 
flavor space. With explicit flavor indices, the Yukawa Lagrangian becomes 


/:Yuk = 

- Y’QR,^UR-yssULUR-MULUR + h.c. ( 2 . 6 ) 
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The mass term of the fermion sector is 




— — (uli U l 




vY' 


0 M + ysU \Ur, 


- vY^jdL^d^J^ + h.c.. 


Hereafter we identify 


Ml, = I 


vY' ' 

0 M + ysU , 


= vY'^ 


where I and J run over 1 to 4. Using bidiagonalisation, the mass matrix transforms as 


U\M^Ur = Diag(m„,TOc,TOt,mT), 

V^M'^Vr = Diag(md,ms,mb), 

through rotations of the quark flavor basis 

Ur = dR = jd^j^, 

where Uj and I?* are 4x4 unitary matrix, and 3x3 unitary matrix, respectively. 

The new CKM matrix i® obtained from the charged current. In weak eigenstates 

Ccc = di lU^ +h.c.. 

Rotating into the mass eigenstates, we get 

Ccc = {'DrY, diw^ + h.c.. 

Uqkm i® defined as a 4 x 3 matrix 

VYKM={ldl)[{VR)Y^U^EV, 
where the explicit form of the new CKM matrix are expressed as 


/I 

0 

0\ 


(Vud 

K. 

Vub\ 

0 

1 

0 

, ^CKM — 

Vcd 

14s 

Vcb 

0 

r\ 

0 

r\ 

1 

r» J 

Vtd 

. T r 

T r 

Vtb 

T r 


Vo 0 0/ \VTd Vts VTbJ 


The unitarity of the transformation matrices U and V implies 

^CKM^CKM = V^E^UU^EV = 13x3, 


(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 


which means that the 3 columns of Vqkm orthonormal to each other. For future use, we can complement 
one column to make up a unitary 4x4 matrix expressed as 


fVud 

Vus 

Vub 

K4\ 

Vcd 

Vcs 

Vcb 

Kc4 

Vtd 

Vts 

Vtb 

144 

\yTd 

Vts 

Vrb 

Vta) 


(2.17) 
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In the scalar sector, the new scalar couples to the SM Higgs boson. The general scalar potential is 






+ 


+ * 


+ T^s 


b^S^. 


(2.18) 


The parameter p-i can be eliminated by a redefinition of the scalar field s —> s + cr. 

The minimization conditions at the vacuum (w, u) are used to eliminate the quadratic coefficients 

1 


“t“ l-^SCp^ “t“ 7 


2 I I 1 \ 2 I 


fig — XsU + fi^U + 7:Xs<pV + 


2 ^ 2u ’ 

The second derivatives of the tree-level potential describe the mass squared matrix of (j) and s: 


mis j 


Diagonalizing the above matrix, we obtain the mass squared eigenvalues 


fisp>v + \s4,uv 

fis<j,V + Xs^jUV 2XsU^ -f fi^u - 


2u 


mh,s 


and the eigenvectors 


1^2 

= jK 


h 

S 


1 

^2 




s) 


COS Lp — sin Lp 
sin ip cos p 


where the mixing angle p is given by 


tan 2p = 




mis - ml^ ■ 


(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 


3 Sources of CP Violation 

In the SM, the CP violation is characterized by the quark-rephasing invariant quantity, the Jarlskog invari¬ 
ant [30] 

Jcp = {ml - ml){ml - ml){ml - ml){ml - ml){ml - ml){ml - ml)A, (3.1) 

where 

H = ImKdKbIC.K; (3.2) 

is twice the area of the unitarity triangle of the CKM matrix. As the SM CKM matrix only has 1 independent 
CP phase, the three unitarity conditions give the same area, which is the only CP violating source. This 
quantity can also be written as [31] 

Jcp = -'pet[Hu,Hdl (3.3) 

where Hu = and Hj, = MdM\ are the building blocks of rephasing invariants. In the picture of 

electroweak baryogenesis, this quantity provides a dimensionless CP violation strength Jop/Tl"^ ^ 10“^°, 
which is too small compared to the typical strength of Baryogenesis rj ~ 10“^°. We observe that both the 
fermion masses and the unitarity triangle A suppress the CP violation. Thus we expect that adding heavy 
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quarks could enhance the CP violating effect via the large fermion mass. In general, this heavy quark can 
be chiral or vector-like. A lot of efforts were performed for both the model of 4th generation quarks [32] and 
vector-like bottom quark [33,34], among which was the study of the enhancement of CP violation. However, 
the simplest 4th generation model was ruled out by the experiment data [35,36]. On the other hand, the 
vector-like quark is still alive and provides enhancement of CP violation in a similar way. 

Let us discuss the CP violation strength in our model. First, let’s count how many independent CP 
phases are there in the model. The unitary condition Eq. 2.16 sets 9 constraints to the CKM elements. 
There are also 6 rephasing redundencies, coming from the 7 involved quark fields modular the total baryon 
phase. Finally, the number of independent matrix elements in the CKM matrix is 

12x2-9-6 = 9, (3.4) 

among which 6 degrees of freedom attribute to real rotations, and the other 3 are independent CP phases. 
The 3 phases can be parameterized as the following 3 rephrasing invariants 

Hi = ImKbF4KT4K4, (3.5) 

B2 = ImVtbV^.VTiVtl, (3.6) 

Ha = linVMVt4V:^, (3.7) 


all of which represents the area of a subtriangle of the unitarity quadrangle formed by the complex numbers 
V^lVt 4 , V*i^Vc 4 , KJj,1/u4 . In the chiral limit where u and c are supposed to be massless, rriu = rric = 0, 
the H’s that involve c and u are not observable, hence we only care about the quantity H 2 . 

Similar to the SM case, those unitarity areas are not the only quantities which appear in the CP violation 
processes. In the case that we are interested in, where CP violation occurs simply via the evolution described 
by the Dirac equation [9], the mass matrices also play the game. Therefore, dimensionful quantities like 
Eq. 3.3 should be used to characterize the CP violation strength. In the gauge basis, the mass terms are 

.^mass — ijUblJ . (3.8) 


We decompose the up-type mass matrix as 



(3.9) 


where M“ and to“ are submatrices with dimension 3x4 and 1x4, respectively. Then we can define 
= M„Mt and as before. We can also define another building block of rephasing invariants 

hu = 

It was investigated that in the vector-like bottom partner model [33], the CP violation is characterized 
by 7 Jarlskog-like invariants (J-invariants). The top partner model should be similar. In the simple case of 
chiral limit m„ = = mg = rric = 0, only one of them is independent: 


J = Im trHdHuKhl 

= — rrij )H2, 


(3.10) 


In this work, we only estimate the CP asymmetry using the Jarlskog-like invariants, and leave more 
detailed study via transport equation for future work. To estimate the strength of CP violation in our 
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model, we need to look at the experimental constraints on the heavy fermion mass and the extended CKM 
matrix elements. Current experiments on flavor physics, such as K and B decay and B — B mixing were 
analyzed in the literature [37, 38] by performing a global fitting on the 4x3 CKM matrix elements using 
68 flavor physics observables [38]. The analysis includes the direct measurements of the CKM elements, 
CP violation in —>■ tttt, branching fraction of the decay K'^ —>■ branching fraction of the decay 

Ki^ —>• Z ^ bb decay, B° — B^ mixing {q = d,s), indirect CP violation in B° —>• J/tpKs and 

Bg —>■ CKM angle 7 , branching ratio of B —>■ Xst^£-~ (£ = e,/i), branching ratio of i? —>■ 

branching ratio oi B ^ constraints from B —>■ K* , branching ratio of B+ —?► 

branching ratio of Bq —>■ {q = s,d), branching ratio of i? —> rP, like-sign dimuon charge asymmetry 

and finally the oblique parameters S and T. The results of the global fitting are shown in Table 5 
and 6 of Ref. [38] for mr = 800 GeV and 1200 GeV. The results suggest that B 2 could be as large as 10“® 
for a TeV top partner More importantly, the enhancement from the heavy top quark mass implies a 
J-invariant of order J < lO^^GeV®. Assuming that the typical energy scale during the baryogenesis is the 
critical temperature of EW phase transition Tf., the dimensionless CP violation strength formulated as J/T® 
needs to be greater than the observed baryon number asymmetry. 

77-10-'°, (3.11) 

which sets an upper bound for the critical temperature Tc < 420TeV. 

Given the possible large CP violation effects, we also need to check the current constraints from the 
non-observation of the electric dipole moment (EDM) of the electron and neutron. The electroweak sector 
of the Standard Model gives an EDM for the neutron of size jd„j ~ 10“®^ — 10“®^ecm. The model with 
an extended quark sector typically gives rise to a quark EDM at the two loop level, thus contributes to the 
neutron EDM. The current experimental limit on the EDM of the neutron is jc?„j < 2.9 x 10“^®ecm (90% 
CL) [39]. It is shown in the literature [40] that the EDM in the model with an extra down-type singlet quark 
and found that the induced neutron EDM is of order 10“^®ecm. Ref. [40] also comments that in the model 
with an extra up-type singlet quark, the down-type quark EDM’s vanish identically at two loop order. Since 
all down-type quarks are light in this case, the leading terms for the Ue EDM is proportional to 
The contributions from the vector-like up-type quark to the neutron EDM are thus even smaller compared 
to the extra down-type quark. 


4 The Scalar Potential 

4.1 A Brief Review of Effective Scalar Potential 

To study the phase transition, we consider the potential of the two scalar fields at finite temperature (see 
Ref. [ 8 ] for review). At the one-loop order, the zero-temperature effective potential in the Landau gauge ^ 

^ Using the results in table 6 of Ref. [38] where the moduli of the V^V elements are estimated by mixing and decay of B 
and K mesons, the Bi quantities can be estimated by Bi < |U^U| x |U^y|. Taking imaginary parts might introduce one or two 

orders of magnitude smaller, but would not ruin the estimation. 

^While the effective potential in the Landau gauge is not gauge invariant, the potential at its minimum is well-defined. For 

concerns about the gauge invariance and a treatment of the gauge invariant effective potential, see Ref. [41]. 



has the form 


Vcw(^, s) = Vtree(^, «) + X/ 


m, (</), s) log —- a 


647r2 




(4.1) 


where rii is the number of degrees of freedom of the particle i running in the loop, with negative sign 
for fermions, and m‘f{(p,s) is the corresponding field-dependent squared mass, defined in Appendix A.3. 
Here Ci are constants that depend on the renormalization scheme, and Q is the renormalization scale. For 
convenience, counter-terms Vct are chosen to preserve the input parameters, like the vacuum expectation 
values (vev) and the masses, from loop; corrections: 

diVcw + K:;t) 


d(j)i 

d^iVcw - 


0i = (0i> 


= 0 , 


Vct) 


d(l)id(l)j 




= m. 


(4.2) 

(4.3) 


where {(f)i) = v,u and mfj are tree level vev and mass squared matrix defined in Sec. 2. This naive treatment 
fails when we consider the Goldstone contribution. The Goldstone boson contribution to the scalar masses 
in Eq. 4.3 is infrared log-divergent due to its zero pole mass. This indicates that renormalizing the scalar 
potential at zero external momenta, as is done in the effective potential calculation, is not a well-defined 
procedure when Goldstone bosons are involved. An alternative on-shell renormalization procedure was 
proposed [15,45,46] to cure the problem, as described in Appendix A.l and A.2 in detail. We extend the 
results of the Ref. [46] to the effective potential with mixture of the Higgs boson and new scalar. Here we 
list the final expression of the zero temperature one-loop effective potential: 

mj{(j),s) 2' 


= ^tree(<^,s) + ^ 


i^G 


647r^ 


Wi (<(<, s) log 


rf(u,zi) 


2mj (u, u)mi {(j), s) 


647r^ 


TOg(<(>,s) In 


mQ(u, u) 




The one-loop thermal corrections to the effective potential at finite temperature T is 


Vthermal(^7F) — ^ ] 27r^ '^B,F 


where 


JB,F{y) = / dx X log 
Jo 


iTe 


rp2 


-\/x^+V 


(4.4) 


(4.5) 


(4.6) 


with the sign — for bosons and -|- for fermions. 

The finite-temperature potential needs to be corrected, due to the infrared divergences, generated by 
bosonic long-range fluctuations called Matsubara zero modes. This can be solved schematically by resumming 
over all diagrams with bubbles attached to the big loop [42-44], which are called the ’’ring diagrams”. This 
leads to a shift of the bosonic field-dependent masses (^, s) to the thermal field-dependent masses 


TOi s, T) = m, ((/), s) + ni(T), 


(4.7) 


where the thermal shifts H^ are defined in Appendix A.3. After resummation, the ring-diagram contribution 
to the effective potential reads 

Vring = - ^ 

2 =B 
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In the SM, this cubic term from the Matsubara zero mode is the only source to induce a thermal barrier 
between a symmetric minimum and a symmetry-broken minimum in the effective potential. It was because 
that all the other terms in SM Higgs sector are quadratic or quartic, which can’t create such degenerate 
minima in one-dimensional scalar space. However, in our model, as shown later, the new dimension in the 
scalar space greatly enriches the possibility, and hence the ring diagram contribution is much less important. 

The total effective potential at finite temperature is the sum of the above terms 

S, T) = S) + Khermal(<^, S, T) + S, T). (4.9) 

For part of the field space, the field-dependent masses of the scalars and the Goldstone bosons can be 
negative, and the non-convexity of the potential would induce an imaginary part that indicates a vacnum 
decay rate per unit volume. However, the real part can still be interpreted as the expectation value of the 
energy density. Therefore we only take the real part of the potential to do the numerical analysis. 


4.2 Approximate Analysis of the Scalar Potential 


In the next section, we perform a numerical study on the full effective potential based on eq. (4.9), and scan 
the parameter space for strong first order phase transition. But to understand the numerical results, we need 
some approximate methods to analyse the complicated potential function. In the high temperature limit, 
the effective potential can be simplified as 


Veff (</>, S, T) ~ Vfcree + ^ 


12tt 


m [m^ 


{4',s,T)] 


3,/2 


- E 


4 


i=B,F 


647r^ 


m, (</), s) In ■ 


I'T^ 


E 

i=B.F 
,2 


24(1 + S,f) 




- 2m, {v, u)m, s) + p{T), 


(4.10) 


where Ug = exp(— 27 ^) for bosons and Op = tt^ exp(— 27 £;) for fermions. p(T) = ^UpT'^ is the Stefan- 

Boltzmann contributions with Up = ns + which is field-independent and can be drop out for our 

purpose. If we series expand the second term in m|((/), s)/T^, the remaining relevant term is another m^T^, 
which, combined with the third term here, gives the main temperature dependence of the full potential. 
The log terms can be absorbed into the running parameters that vary little within the energy scope of our 
discussion. Now the effective potential is simplified as 


VMs,T) ~ 


me. 


1/2 




Sir 




24(1 + ,5,f) 


m^(</.,s)r2+p(r), 


(4.11) 


where Vq is the zero-temperature part of the potential, being the tree-level potential plus loop suppressed 

1/2 

corrections. Taking the approximation cd /Stt « 1/24 for all bosons, we obtain the temperature-dependent 
terms 


Heff (/>, S, T) ~ Ho + [Hs + + Cs^] , 


where the coefficients are 
A 
B 


1 

12 

1 

12 


[Aps4> + 2/i3 -I- 6 ysM], 

( 6 A 0 -I- As0/2) -I- ^(3g^ + 5 ^) + + v'^) 


(4.12) 


(4.13) 

(4.14) 
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c 


Y2 [(3^s + 2As0) 3y^] . 


(4.15) 


Note that the fermion Yukawa couplings are also contributed to the configuration of the scalar potential. 
This simplification results in a polynomial form of the potential, which is convenient for our analysis. There 
are sometimes significant errors for this simplification, but we will show that the qualitative analyses based 
on this polynomial potential explain many key features of the numerical results we obtained. 

We will analyse the scalar potential at the moment of the phase transition. More specifically, we are only 
interested in the stable vacuum, i.e. the global minima of the potential, which is degenerate at the critical 
temperature of phase transition. They are defined as (0 ,Ms) and {vc,Ub), representing the symmetric vauum 
and symmetry breaking vacuum respectively. In the following, we introduce two kinds of method to estimate 
the properties of these vacuum configurations. 


4.2.1 Barrier Width Estimation 


We used to express the potential in the cartetian coordinates of the field space {(j), s). For the analysis of the 
phase transition, it is also convenient to utilize the polar coordinates {p,a) [18,22,47] 


P = + (s- u) 2 , COSO = 

y'p + (S - W) 

from some shifted center (0 ,m). When u = Us and 

Vc 

cos On = — , 

y/v'^ + (Ub - UsY 


(4.16) 


(4.17) 


The potential has degenerate minima along the p axis, p = 0 being the symmetric one, and p = p being the 
symmetry breaking one. Hereafter we will use the notation Ca = cosoq and Sa = sinoo for short. Similar 
to the SM, we can employ the following form of parameterization of potential 

Heff (p, T) \d{T^ - T^)p^ + £p3 + ^p\ (4.18) 

where the coefficients are constant for a simplified polinomial potential. The coefficients D, £ and A are 
functions of the model parameters and Us, Sa,Ca. We impose the following condition 


Ves{p,Tc) = \p^{p- pf, 


and obtain the non-zero vacuum value 


2 £ 



If we neglect the zero-temperature loop corrections, they can be expressed as 


£ — [(Msc/) T As0'Us)Cq,/ 2 -f (/r3/3 -f As'as)sQ,j , 

A = A^c^ -l- Aj;S^ -l- Ajs^c^s^. 


(4.19) 


(4.20) 


(4.21) 

(4.22) 


The p determines the width of the barrier in the scalar potential, which is what we finally concern. But this 
analysis does not give estimations of Us,Sa,Ca, which brings us to the next tool. 
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4.2.2 Stationary Point Search 


Although the shape of scalar potential has been studied in Ref. [19], we provide a detailed and systematical 
recipe to describe the shape of the 2-dim potential. We summarize our results in Table 1, which could be 
used to understand the numerical studies of the phase transition in the next section. 

First, let’s write the potential 4.12 as the following 


Vee{(t),s,T) 


^ 1 
4 ^2 



+ jj-sipS 



X 


~ - 2 

.. 4 I 3 /^s 2 I ~3 


(4.23) 

(4.24) 


where all the tilded couplings are supposed to depend on temperature, logarithmically or quadratically. Since 
the former are negligible within the energy scope that we are interested in, we only need to focus on the 
quadratic temperature dependencies 



= AT^, 

(4.25) 



(4.26) 

As 


(4.27) 


while the other couplings are mainly their zero temperature values. We learn that a temperature around 
100 ^ 150 GeV is usually smaller than the other massive parameters, therefore even the quadratic temper¬ 
ature dependencies are still insignificant at this range of temperature. As a result, the coefficients roughly 
satisfy 


Ixl < lAsI ~ ImsI- 


(4.28) 


Neglecting the linear term further implies a useful corollary, that there is always a stationary point sitting 
around the original (0, 0), even until the temperature reaches 100 ~ 150 GeV, the typical values of the critical 
temperature in our model. Only at higher temperature when y becomes important, will this stationary point 
gradually move away. This corollary is verified by the parameter scan, and is essential for the explanation 
of some of the transition patterns. 

Now let’s do a thorough search of the possible stationary points in the potential. First we notice that 
the condition for extrema consist of the following two curves: 

ocp os 


and the vacuum must be at one of the intersections between the two curves. Let us call them the (j) curve 
and the s curve respectively. We will describe the shape of these two curves, and try to find some rules that 
the possible degenerate vacuum points should follow. 

The (j) curve consists of a trivial line (f) = Q and a quadratic curve 


, r/2 

\ _ "2 I ^S(p 


+ =~4 + 


^S(p 


2A 




(4.30) 


Obviously, symmetry-broken minimum, if there is any, must be on the quadratic curve. 
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For the case of As 0 > 0, the quadratic curve is an ellipse, centered at point (0, s* = —with 
size decreasing as —T^ due to the term. At some high temperature, the ellipse shrinks to zero, and no 
non-zero vacuum is allowed, hence the symmetry must be restored. If Ag^ < 0, the quadratic curve is a 
hyperbola also centered at (0, s*). At high temperature, the curve will move away to infinity. As long as the 
potential is still bounded from below, the minimum cannot be on the hyperbola, so the symmetry must be 
restored. 

The equation for s curve is a cubic polynomial of s 

Ags^ + Ass" - [nl - ^Ag^(^2^ s+(^f + = 0 (4.31) 

The discriminant of this polynomial turns out to be a polynomial of the </> field: 

A((()) = a(j)^ + + c(t>^ -\- d (4.32) 


where the coefficients are 

® ~ 2 

u ^ \2 ~2 I o\ \2 ~2 I 9r r _ _ 27~2~2 

^ ^ '^sdsip! 

c = Ag^/ig/ig 6 AgAg,^/ig 9 Ag/i 3 /ig/ig,^ 

d — (fil + 4Ag/ig)/ig. 


(4.33) 


The insignificant x has been taken to be zero as an approximation. This polynomial shows the number of 
points on the s curve at a specific (j) value: if the polynomial is positive, then there are three points at this 
(j); if the polynomial vanishes, two of the points are degenerate; if the polynomial is negative, there is only 
one point at this (j). Therefore, the variation of the sign of this discriminant with respect to (j) could give the 
key features of the shape of the s curve. A systematic classification is provided in Table 1. 

In the table, the middle column shows how the sign of A(^) varies with (j) (red indicates a possitive 
discriminant while blue indicates a negative one). The typical shape of the corresponding s curve is given 
in the last column. The first colume shows the condition for these types - when the A' = — 4ac^ — 

27a‘^d^ + ISabcd is positive, all the signs of a,b,c,d should be specified; when A' < 0, only a,d need to be 
specified. 

The upper half of the types (A,B,C,D) have positive coefficient a, and 3 asymptotic values at large (j). For 
the lower half (E,F,G,H), where a < 0, there is only one asymptotic value at large 4>. Also note that the left 
half of the types (A,C,E,G) have positive coefficient d, which leads to multiple branches (a branch is defined 
by a topologically connected part of the curve) in the small (j) region. The right half, however, have negative 
d and only one such branch. We name the branches present at small (j) as ” relevant branches”. Usually, the 
’’irrelevant branches”, which only appear in large </> region like in type (B,C,D), are not important in our 
analysis. This reduces the set of types by identification B ^ H, C E and D ^ F. 

From the point of view of Lagrangian parameters, we find that the most important coefficients a and 
d are controlled by Ag^ and -I- 4Ag/Xg. The latter changes with temperature mainly through Eq. 4.25. 
Moreover, the asymptotic value of the s curve at large (j) (the middle one when a > 0) coincides with the 
center of </> curve s*. 

In the next section, we will show how the analysis of the 4> curve and s curve could help classify the phase 
transition patterns. 
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Table 1: The s curve shows the first derivative of the scalar potential along the s direction ^ = 0. A' is 
defined to be the discriminant of A((/)) as a cubic polynomial of (j)^. It is positive unless specified. 


'I’ype 

{a,b,c,d) discriminant A ((;A) 


A 

(+i +1 +i +) 


C 

(+i +i +) 

(+i +i +) 
(+i “I +) 


E 

(—, +1 +, +) 
(“I +i +) 

(“I “I +) 



'i'ype 

(a, b, c, d) 


B 

( + 1 +: +: “) 
( + i +: “) 

( + i “) 



discriminant A((/)) 



*[GeVl 


s curve ^ = o 



*[GeV] 


G 

(“> +. “I +) 



*lGeVJ 


'I'ype 

(a, d) discriminant A(0) 


A 

(+i +) 

A' < 0 


E 

(“I +) 

A' < 0 


s curve 


'I'ype 

(a,d) 


av 

ds 


B 

(+i “) 

A' < 0 


H 

A' < 0 
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5 Strong First Order Phase Transition 


Once we obtain the full effective potential, we could investigate how the vacuum state evolves with temper¬ 
ature. At each temperature, we find the true vacuum by looking for the global minimum of the potential in 
the (0, s) field space. It is known that at zero temperature, the global minimum of the scalar potential is at 
(uo, uq), with vq = 246 GeV and uq as an input parameter. After we turn on and increase the temperature, 
we track the position of the global minimum, seeking the sign of a phase transition. In the context of effective 
potential, the fields are defined to be the first derivative of the free energy with respect to the corresponding 
particle source, and thus acts like an order parameter of the phase transition. If the transition of the field 
values between the two phases is continuous, it is a second order phase transition. Otherwise, if there is 
a discontinuous variation of the fields, it is a first order phase transition. The first order phase transition 
proceeds by bubble nucleation of the broken phase at around the critical temperature. The bubbles grow 
and coalesce, and finally turn the whole universe into the broken phase. 

As discussed in the introduction, to have successful baryogenesis, it is essential to have a strong enough 
first order phase transition, so that i?sph(T'c)/T'c > 45 inside the bubbles. It has been shown in literature [17] 
that for a singlet extended model the sphaleron energy is approximately proportional to the (j) vev: illsph(7c) ~ 
Vc- Moreoever, the bubble expansion and wall velocity in the singlet extended model have been discussed in 
Ref. [48]. In our model there is a similar scalar sector. Thus the discussions about the bubble expansion and 
sphaleron process in literatures are also applied to our model. Therefore, similar to singlet extended model, 
we add the following criterion to our scan to pick out the events of successful baryogenesis: 

C=g>l. (5.1) 

We also perform a consistent check by calculating the sphaleron profiles and the sphaleron energy at the 
critical temperature numerically. 

To determine the parameter region in which the strong first order phase transition could happen, we 
perform a random scan over the parameter space. The procedure is the following. We have quite a few 
independent parameters 


Af, y , (5.2) 

in addition to the singlet vev uq in the zero temperature. The parameters X^f, and yt are determined by 
the Higgs mass and the top quark mass that are already known. We choose the input parameters from the 
ranges ^ 


< 1.5, 0 < A, < 2, \y'\ < 1.5, |y«| < 1.5, (5.3) 

|wo| < 600, l/ral < 800, \fis,p\ < 1000, 0 < M < 1200. (5.4) 

Given the input parameters, the full effective potential in the field space is calculated. Then for each 
temperature, we utilize the MINUIT subroutine [49] to find the global minimum of the effective potential. 
As the temperature increases, we track the change of the global minimum at each step in our numerical 

^ When the singlet vev uq was chosen to be larger than 600GeV, the data with critical temperature larger than 200GeV, and 
a Vc even larger, accumulates. They are mostly the case IIIB as introduced later, and are not our focus in this paper. That’s 
why we chose a smaller range for uq. 
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scan. Additional care should be taken: if the minimum moves to very large field values, it may indicates the 
vacuum instability. If the global minimum becomes the symmetric one (0, (s)) at certain temperature, we 
perform a fine scan near the temperature until we hnd the critical temperature Tc and the corresponding 
vevs {vc,Ub) in the broken phase. After obtaining the Tc and Vc, we use the washout condition to pick 
out the strong first order phase transition, eliminating the data points that has C < 1- We randomly scan 
10® parameter points, among which 25818 parameter points pass all the requirements. Figure 1 shows 
the distribution of the successive data points in the Tc — Vc plane. From the figure we notice that for the 
parameter region we scanned, the critical temperature is typically less than 200 GeV, while Vc is smaller than 
its zero-temperature value vq = 246 GeV. As expected, the higher the critical temperature is, the smaller 
the (j) vev gets to before the transition, and this correlation is clearly shown in the figure. Furthermore, this 
range of critical temperature is quite safe from the bound of the GP violation strength that we discussed in 
Sec. 3. 



Figure 1: The allowed value of the critical temperature Tc versus the (j) vev Vc at the critical temperature from a 
random scan over the parameter space. The scatter points are selected to satisfy the SFOPT in our random numerical 
scan. The color palette on the right shows the density of the scatter points in one GeV interval. 


5.1 Phase Transition Patterns 

The patterns of phase transition are described by the critical values of the scalar fields during the transition, 
dehned as (0, Us) and (uc, Ub) in the previous section. In the left panel of Figure 2, we show the correlation 
between Ub and Us, indicating the s value jump during the phase transition. We notice that there are several 
distinct allowed regions: Region I: the vev Us is very close to zero while Ub is non-zero; Region II: the vev Ug 
is non-zero linearly correlated with the Ub] Region III: both the vev Ug and Ub are scattered over the second 
and fourth quadrants. The Region III can be further classihed by the correlation between Ug and uq, which 
characterizes the evolution of s value before the transition. This is shown in the right panel of the Figure 2. 
Typically, the vev uq and Ub are strongly correlated, because there is not much time for s to move very far 
without jumping. However, there are also cases where Uq and Ub are quite uncorrelated. This indicates that 
between the zero temperature and the critical temperature, there should be another phase transition, though 
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Figure 2: The allowed value of the contour {ua, ui,) (left) and {ub, uq) (right) from a random scan over the parameter 
space. The color palette on the right shows the density of the scatter points in one GeV interval. 


the Z 2 symmetry is not restored yet. Further analysis shows that the region that the uq and Ub are not 
linearly correlated only appear in Region III. So we separate the Region III into two subregions: one-step 
phase transition (Ilia) and multi-step phase transition (Illb). Hopefully, these four regions of parameters 
correspond to four different phase transition patterns. 



Figure 3: The allowed value of the contour {£, X) (left), (he, Vc) (middle) and (£, ps<#> + Xa<i>Us) (right) from a random 
scan over the parameter space. The color palette on the right shows the density of the scatter points in one GeV 
interval. 


In order to help estimating the parameter preferences, we also utilize the polar coordinates of the scalar 
fields as in Eq. 4.18. In terms of the non-zero vacuum value p, we obtain the (f> value for the approximate 
potential: 

- - 

Vc = pCa = -^Cq. (5.5) 

A 

A necessary condition of having strong first order phase transition is that is positive and large. However, 
the expressions of the £ and A cannot be expressed as functions of only the model input parameters. They 
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also involve the information of the features of the phase transition, such as the angle a and the symmetric 
minimum Ug. Therefore, we would like to further utilize the classification of the transition patterns to get 
more information on these features. 

Several approximations can be made here, according to the correlations between the parameters that we 
found. First, in the Region II and perhaps also Region I, where Ug and Ub differ only a little, we thus assume 
that a <C 1. For Region I, Ug is very small in comparison with all the other massive parameters. After these 
approximations, the relations between Vc and the couplings will be more manifest. In Figure 3 we plot the 
comparison between the Vg and the true Vc that we obtain from the scan (left), and the correlation between 
£ and + Ag^Ug (right). The left panel tells us that in most cases Vc is a good approximation to the true 
value. The right panel shows the correlation between £ and its main part under the assumption of a <C 1. 

In order to analyze the four regions in detail, we also have to rely on the shape of the scalar potential. 
The (j) curve and s curve at the critical temperature help us understand the origins of the different patterns, 
as they determine the distribution of the potential minima. The s curve could have one, two or three relevant 
branches, depending on the types as shown in Table I. Therefore, it could have one or three intersections 
with the line (p = 0. We call the minima along the line p = 0 ’’symmetric minima”, while others are 
called ’’broken minima”. Next, we notice that the broken minimum must be the intersection between the 
quadratic branch of the p curve and one of the s curve branches. This s branch must have an intersection 
with the s axis, which may be a symmetric minimum. We would like to define the barrier between the 
broken minimum and this symmetric minimum as a ’’single-branch barrier”, while those between minima 
on different branches are called ’’inter-branch barrier”. We discovered that the single-branch barrier usually 
has much smaller width along s direction than the inter-branch barriers, due to the limited stretch of the s 
curve along the s direction. It implies that a transition through the single-branch barrier should have closely 
related Ug and rtf,. Finally, in terms of the above features of the shape of potential, let us discuss the four 
phase transition patterns in detail. 

5.1.1 Pattern I: Single-branch barrier transition, with Wg ~ 0 

According to the discussion in Sec. 4, there are three cases for the small </> behavior of the s curve, one of 
which is usually negligible. In the other two significant cases, the one with three roots will be discussed in 
the next part. Let’s focus on the other case, which has only one relevant branch that intersects p = 0 at 
around (0, 0). It indicates that this case mainly corresponds to the Region I of the parameter scan. 

In Figure 4, we show the allowed values of the broken minimum (vc, Ub) in this case, while the symmetric 
minimum is always at (0,Mg ^ 0). The value of Vc is upper-bounded by its zero-temperature value vq = 246 
GeV, implying a decrease of p before the transition, and is also lower-bounded at about 100 GeV by the 
condition ^ > 1. The distribution of the broken minimum clearly sketches the shape of the s curve. 

Let’s see what can be inferred for the model parameters. Without loss of generality, we choose the 
benchmark points with only positive uq, to investigate the sign preferences of other parameters. Here are 
the observations: 

• As shown above, the s curve is a single branch curve across the point (0,0), which indicates that the 
zero-temperature vev ug, which is on the same branch, should also be small, but not zero due to the 
bend of the curve. 
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Figure 4: The allowed values of the broken minimum {vc, ui,) in pattern I. Here the phase transition happens between 
the symmetric minimum (0, 0) and the broken minimum. The color palette on the right shows the density of the 
scatter points in one GeV interval. 


• In the light of the previous discussions, there are several approximations we can employ for £: Ug « 0 
and a <C 1. As a result, the only important term in £ is the fj.s<p term: 

£ ~ (^■^) 
In order to get a large and negative £, large and negative ^s(p should be favored. 

• We would like to argue that positive Ag^ leads to shapes of the curves that are much more favored by 
the strong first order phase transition. One may notice that for negative both </> curve and s curve 
are hyperbola-like. Two hyperbolas could not make the twisting intersection needed for the existence 
of degenerate minima. Although the parameter flg^ causes a deviation from perfect hyperbola of the s 
curve, it’s still harder for such case to have first order phase transition. Therefore, positive Ag^ should 
be strongly preferred. 

• The one-branch condition for s curve requires that the coefficient d in the polynomial A(^) is negative, 

-2 

thus jll < — . Large would compress the parameter region that satisfies this condition. Thus in 

this pattern, we expect that small /i§ is favored. 

Fortunately, our numerical results from the scan do exhibit the above features, as shown in Figure 5. 

Figure 6 and Figure 7 show the variations of the (f> curve and s curve with temperature. In both 
figure, the first diagram represents the configuration at zero temperature. The second, third and fourth 
diagrams represent the configuration below, at and above the critical temperature. In the diagram at critical 
temperature, an arrow was drawn to show how the phase transition happened, in the point of view of 
increasing temperature. Similar figures will be given for other patterns of phase transition later. 

The difference between the two figures is whether the barrier exists or not at zero temperature. In 
Figure 6, only broken minimum exists at zero temperature, then the symmetric one is developed during the 
heating. In Figure 7, the symmetric minimum already exists at zero temperature, but with higher potential 
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Figure 5: In the pattern I with uo > 0, the allowed values of the model parameters uo, and ^3 are shown. 

The color palette on the right shows the density of the scatter points in one GeV interval. 


than the broken one. The existence of the barrier at the zero temperature is a new feature for the singlet- 
assisted electroweak phase transition, which does not happen in the traditional electroweak phase transition 
where the barrier must be thermally induced. This could be attributed to the existence of the term /is^, as 
shown in the expression of £ where only the is important in pattern I. In other patterns, other terms 
in £ could also contribute to the barrier. Unlike the SM, these contributions to £ don’t require a non-zero 
temperature, hence zero-temperature barrier can exist. 

5.1.2 Pattern II: Single-branch barrier transition, with Us separated from 0 

Suppose the s curve has at least two relevant branches. Assuming /us is small, we only consider the case when 
/Tg > 0. As mentioned earlier, there are two symmetric minima, one positive and the other negative, while 
the stationary point around (0,0) must be a saddle point. Let’s focus on the single-branch transition for 
now, which means that the symmetric minimum (0,M)f*’®'') on a different branch from the broken minimum 
has higher potential than (0, on the same branch as the broken minimum. The other case will be 

discussed later. 

Now that the symmetric minimum chosen by the phase transition is other than (0,0), we are convinced 
that Us = can’t be close to 0. That being said, the gap in Figure 8 is well understood. Meanwhile, Us 

is still strongly correlated with Ub because they are on the same branch. It is precisely the characteristic of 
the transition pattern II. 

Let’s investigate the parameter preferences in this pattern. As the whole branch of s curve under consid- 
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Figure 6: In the phase transition pattern I, the s curve (dashed blue) and (p curve (solid green) in the different 
temperatures: zero temperature, below critical temperature, at the critical temperature, and above critical tempera¬ 
ture. A thermal barrier is developed during phase transition. The color palette on the right shows the density of the 
scatter points in one GeV interval. 






Figure 7: In the phase transition pattern I, same as the Figure 6, but a tree-level barrier is developed. 


eration is apart from s = 0, it is natural to have uq also apart from 0. Actually, all of the uq, Ub and Ug are 
similar, which means that s field is seemingly irrelevant with the phase transition, although it contributes 
a crucial coupling that plays important role. The preference of the is similar as pattern I. A big 
difference from pattern I is that positive As^ is strongly suppressed. One can understand it by looking at 
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Figure 8: The allowed values of the broken minimum (nc, Ut) in pattern II. The color palette on the right shows the 
density of the scatter points in one GeV interval. 






Figure 9: In the pattern II with uo > 0, the allowed values of the model parameters uo, ^s(i> and /rs are shown. 
The color palette on the right shows the density of the scatter points in one GeV interval. 


£, which in this case has two significant terms: and Ag^Ug, the latter becoming significant due to the 

non-zero Ug. As we need a large negative £, negative Ag^ is preferred. Therefore, the advantage of positive 
Ag 0 is much weakened in pattern II. Another interesting feature is that is mostly negative. This is essential 
to guarantee that among the two symmetric minima, the one with positive s value has a lower potential. 
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and tends to be the one chosen as the vacuum state at high temperature. Otherwise, inter-branch transition 
should happen. 



s[GeV] s[GeV] 


Figure 10: In the phase transition pattern II, the curves are same as the Figure 6. 


The procedure of the phase transition is outlined in Figure 10. The only difference from pattern I is 
the existence of the symmetric minimum on the other branch. This minimum must already exist at zero 
temperature in this pattern, which induces a barrier that is impossible in traditional electroweak phase 
transition. As for the single-branch barrier, similar to the pattern I, it can be thermally induced, or already 
present at zero temperature. 

5.1.3 Pattern Ilia: Inter-branch barrier transition 

The setup in this case is similar to the pattern II, except that the phase transition occurs across the inter¬ 
branch barrier. 

As discussed in pattern II, /ra determines the relative height between the potential at the two symmetric 
stationary points (0,u|“®*®) and (0 ,m™*^®’'), where u|™®'® and are the scalar s vev in single branch and 

inter-branch. It tends to be negative when we want a single-branch transition towards (0, u|™®^®), and for 
the same reason, it prefers positive values when we require an inter-branch transition towards (0 ,m“*®'’). 
However, if /ia is too large, it would be harder for us to get a large and negative £, as in this case Sa is 
not small, and hence the /ias^ term in E is not suppressed any more. In sum, we should have a small and 
positive ^a in pattern Ilia, which is justified by the Figure II. 

The phase transition is shown in Figure 12, by exhibiting the variations of the (j) curve and s curve with 
temperature. The only difference from pattern II is that the phase transition happens across the inter- 
branche barrier. Nevertheless, as the (0, need not be a minimum, no twisting intersection is required 

for the two curves, and hence there is no strong preference for Ag^. 
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Figure 11: In the case Ilia with uo > 0, the allowed values of the model parameters uo, ^ 34 , and (13 are shown. 
The color palette on the right shows the density of the scatter points in one GeV interval. 






Figure 12: In the phase transition pattern Ilia, the curves are same as the Figure 6. 
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5.1.4 Pattern Illb: multi-step transitions 


In our investigation of the phase transition, we focused on looking for the critical temperature when the (f) vev 
jumps from a non-zero value to 0. We didn’t investigate whether there was a jump below this temperature. 
Nevertheless, we learn from the correlation between the zero temperature vev uq and the critical temperature 
vev rtf,: if there was any big difference between them, we would expect another phase transition. In our scan, 
we find a rough linear relation between uq and Ug (shown in Figure 13), implying that the zero temperature 
vacuum and the above-critical temperature vacuum are on the same branch. 






Figure 13: In the case Illb with uo > 0, the allowed values of the model parameters (uo, Ub), Vc, and ps are 
shown. The color palette on the right shows the density of the scatter points in one GeV interval. 

There are many possible ways of realizing this. One of them would be that a elliptical (j) curve is placed 
across multiple branches of the s curve, forming several broken minima, as shown in Figure 14. The global 
minimum may transit from one branch to another, and jump back later, causing the UQ,Ub,Us correlation 
previously described. 

Unlike other cases, one could find more events with small values of in this pattern, indicating a small 
s* and thus a cj) curve centered near ( 0 , 0 ), which may be more likely to intersect with multiple branches 
of the s curve. The signs of the parameters can be inferred from the comparison of potentials among the 
three points: negative /is is preferred for a final minimum at (0, Ug >0), and positive is preferred for a 
minimum at (vc,Ub < 0). Figure 13 precisely shows these features. 

As the Figure 1 shows, events concentrate at the region Tc < Vc < 246, with Tc typically around lOOGeV. 
In this pattern Illb, however, we found that the critical temperature tends to be large, and Vc is even larger, 
which is very different from the other patterns. Large critical temperature leaves sufficient space for a second 
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Figure 14: In the phase transition pattern Illb, the curves are same as the Figure 6. 


phase transition below it, but may be harmful for a strong enough CP violation strength because it appears 
in the denominator of the strength with 8 powers. Our estimation of the upper bound for Tc pretty much 
rule out this pattern as candidate of EW baryogenesis. Another interesting feature is that this pattern 
becomes more and more likely as Uq becomes larger. Test scans showed that a larger range of uq would 
lead to dramatically more events with pattern Illb. In order to evade the interference from these ’’complex” 
situations, we chose a smaller range of uq for the scan. 

In summary, the above four cases are the typical patterns in our results of the parameter scan. The first 
two cases are the most concerned when uq is within the typical range of energy scale for new physics, in that 
they account for 90% of the data from our scan. The pattern Ilia is overall very rare, but the pattern Illb 
would become favored at large uq. 

5.2 Numerical Results 

Given the analysis about the transition patterns, we would like to present our numerical results on physical 
parameters. Since we performed a random scan on the parameter space with flat prior, the distributions of 
the parameters should reflect the preference in the model to have the strong first order phase transition. 

Figure 15 shows the allowed values of the broken minimum (vc,Ub) and all eight input parameters. This 
is basically a summary of the detailed discussion in the previous subsection. To complete our discussion. 
Figure 15 also shows the parameters in the fermion sector. We learned that the fermion sector is almost 
irrelavent to the strong first order phase transition. 

We also show the allowed region for the derived parameters in Figure 16. We found that the value of the 
is bounded from a minimum value of about 0.1, and peaks at around 0.5. This indicates that a Higgs 
self-coupling stronger than that in the SM is expected in the singlet scalar extended model. The as 
a controller of the intersection between the s curve and the s axis, has been discussed in the last section. 
Obviously, for patterns I and II, it tends to be negative and positive respectively. Interestingly, Figure 16 
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Figure 15: The allowed values of the broken minimum (vc, Ub), and all eight input parameters allowed by the strongly 
first order phase transition. 


shows that it has no specific sign preference over all. The parameter characterize the intersection between 
the (j) curve and the (j) axis. There are non-zero intersecions only when /i^ > 0. The high peak in Figure 16 
represents the pattern I in which (/) curve roughly go across the origin (0,0). The large tail on the negative 
side represents pattern II, in which the (j) curve is usually far away from the (j) axis. 

Besides the input parameters in the model, we also obtain the favored region of the physical observables, 
like the masses and mixing angles of the new particles, in light of the strong first order phase transition. 
In Figure 17, we present the two dimensional contours of the physical parameters (ms, sin a) in the scalar 
sector, and those of {rriT, sin0) and (M, mr) in the fermion sector. It is shown that the scalar with its mass 
around 500 - 1000 GeV and a medium mixing angle is favored. We recognized the feature that small mixing 
angle are disfavored, as expected from the fact that the scalar needs to couple with the Higgs boson to render 
strong first order phase transition. This favored region is compatible with that allowed by vacuum stability 
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Figure 16: The allowed values of the derived parameters /ij and in GeV. 



Figure 17: The allowed values of the contours of the physical mass and couplings (m5,sin(/5) in the scalar sector, 
and (mT,sin0i) and (M, my) in the fermion sector. 


criteria [23]. Unlike the scalar mixing angle, the fermionic mixing angle can be very small, which indicates 
the decoupling between the new fermion and the phase transition criteria. We expect that the precision data 
and the Higgs data will put stronger constraints on the fermion sector. Finally, the strong correlation in 
the last panel implies that the mass of the new fermion is mainly controlled by the Dirac mass term, thus 
M + ysUQ « M, which may result from a relatively small region of uq that we chose. 


6 Implications at the LHC 

From the numerical results, we found that different transition patterns exist, among which parameter pref¬ 
erence are different. Regarding to the physical parameters, we note that a new scalar boson with 500 ^ 800 
GeV mass and medium mixing angle are favored. 

In this section, we check whether the strongly FOPT parameter region is still allowed by the current 
experimental constraints, such as the oblique corrections S, T, Higgs coupling measurements, and direct 
collider searches. 

The oblique corrections S', T put the tightest constraints on the fermion mass and mixing angle. For a 
singlet vector-like fermion, the Zhh measurement is less stringent than the oblique correction T. Here we 
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collect the results [23] on the oblique parameters into boson-loop contributions T 5 , Ss and fermion-loop 
contributions Tp, Sp. For the fermion-loop contributions, the NP effect is only involved in the vacuum 
polarization amplitudes where the top quark and heavy top quark are in the loop. Subtracting the SM 
contributions due to the third generation quarks 
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we arrive at the expressions from the fermion contributions 
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Similarly the contributions from the Higgs and scalar loop are 
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We note that the scalar contributions are much smaller than the fermion contributions. Thus the constraints 
from S, T on the scalar mass and mixing angle are quite weak. 

On the other hand, the Higgs coupling measurements at the LHC put the tightest constraints on the 
scalar mixing angle. In our model, due to mixing between the Higgs boson and the heavy scalar, all the 
tree-level Higgs couplings are modified as 


NP SM NP SM 

9hff = cos ipglff, ghvv = cos ^glvv ■ 


(6.9) 


At the same time, the loop-induced Higgs couplings to the photon and the gluon are also modified by the 
new contribution from the vector-like fermion loop. So the Higgs couplings to the photon and the gluon are 
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where Tp = and ghtt,9hTT are the Higgs couplings to the fermions, given in Ref. [23]. In the Higgs to 
photon and gluon process, there is almost no fermion mass dependence if the heavy vector-like fermion is 


29 
































much heavier than the top quark. Therefore, the Higgs coupling measurements put constraints on the scalar 
and fermion mixing angles. 

The direct searches on the new fermion and new scalar boson at the LHC also put tight constraints 
on their masses and couplings. At the LHC, the vector-like quark can be produced in pair through QCD 
production pp TT, or be singly produced via electroweak process pp —>■ Tb. In our paper, the up-type 
vector-like quark predominantly couple to the third-generation quarks through T —>■ tZ, th^ Wb. From an 
updated CMS analysis [50] which uses the 8 TeV data collected up to integrated luminosity of 19.5 fb“^, 
the lower limits on the mass rriT are set to be around 687 — 782 GeV. 

The heavy scalar S is CP-even and has the same quantum number as the SM Higgs boson. The search 
limits on the high mass Higgs boson at the LHC could be recasted to put constraints on the mass and 
the coupling of the heavy scalar. The production mechanism is quite similar to the SM Higgs boson. The 
dominant channel is the gluon fusion channel gg —>■ S. The decay channels of the S boson include 

S WW, ZZ, S^hh, S ^ ti, (6.12) 

and S ^ tT only if the heavy top is much lighter than the scalar S. Other decay channels, such as S' —)■ 77 / 53 , 
S —>■ //, where / is the fermion other than the top quark, are negligible. Using the updated analysis [51] 
in the S —)■ WW and S —>■ ZZ decay channels from CMS, the mass range between 145 GeV and 1000 GeV 
has been investigated. According to the analysis, if the high mass Higgs boson has the same coupling as the 
SM, the mass range between 145 GeV and 710 GeV are excluded at the 95% CL. In our model, we recast 
the exclusion limits to the constraints on the scalar mass and mixing angles. 
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Figure 18: The allowed parameter contour (mg, simp) (left) and (mT,sin0L) (right) in light of the strong hrst order 
phase transition. The constraints from the S, T parameters, Higgs coupling measurements, and direct LHC searches 
are shown as the exclusion lines. 


The numerical results on the constraints from the S, T parameters, Higgs coupling measurements, and 
direct LHC searches are shown in Figure 18. From the Figure 18 (left), the parameter region with the mixing 
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angle sin a greater than 0.43 has been ruled out by the current Higgs coupling measurements. The direct 
LHC searches also exclude the scalar boson with mass less than 500 GeV. We expect that the future Higgs 
coupling data put more stringent constraints on the mixing angle sin(/j, and thus put stronger limit on the 
favored region by the strong first order phase transition. On the right panel of the Figure 18, we note that 
the S, T parameters can only exclude very small region which is favored by the strong first order phase 
transition. The direct LHC searches could exclude the top partner with mass less than 700 GeV. There are 
large available parameter regions in the fermion sector. Therefore, the CP violation rate from the fermion 
sector is adequate to generate the needed baryon number asymmetry. 

Finally, we expect that the future Higgs data could explore the parameter region on the {ms, sin ip) 
contour. Furthermore, from the Figures 15 and 16, we note that the large scalar coupling and moderate 
Xfj) is strongly favored. We should be able to explore the scalar trilinear coupling Ag^ and A^ at the high 
luminosity LHC. If the trilinear couplings are enhanced compared to the SM Higgs self-coupling, the Higgs 
pair production cross section should be larger than the SM value. Through the Higgs pair production 
process pp ^ h/S —>■ hh, we could extract out the trilinear couplings from the production cross section 
measurements. 


7 Conclusions 

We investigated the necessary conditions to realize the electroweak baryogenesis in a scalar-assisted vector¬ 
like fermion model. In the fermion sector, the extended CKM matrix provides additional sources of the CP 
violation effects, parametrized by Jarlskog-like invariant. We found that the CP violation rate is greatly 
enhanced by the heavy mass of the new fermion. With the flavor constraints on the extended 4x3 CKM 
matrix considered, we estimated the CP violation strength, which turns out to be adequate for the baryon 
number asymmetry. 

We focused on the one-loop, finite-temperature effective potential in our model and its implications on 
the electroweak phase transition. Unlike the case of the SM, the new scalar extends the field space in which 
the phase transition occurs. In the two-dimensional field space, we have more possible ways of constructing 
barriers between minima. We utilized the shape of the derivative of the potential: s curve and the cj) curve, 
as a tool to analyse the two-dimensional effective potential. 

The first order phase transition occurs in the form of bubble nucleation of the symmetry broken phase. 
The sphaleron decoupling criteria ^ > 1 is used in this model, to prevent the baryon asymmetry 

generated in the symmetry broken phase from being washed out by sphalerons. We performed a parameter 
scan over the 8 independent model parameters, and obtained the allowed parameter region which could have 
strong first order phase transition. 

According to the different regions in the {ub, Ug) contour at the critical temperature, transition patterns 
are classified into four patterns: single-branch barrier transition (pattern I or H, with or without the existence 
of multiple relevant branches in the s curve), inter-branch barrier transition (pattern Hla) and multi-step 
transition (pattern HIb). For the single-branch barrier transition, the large trilinear mass term /ig<^ is favored 
because the width of the barrier is strongly related to it. However, small /ig^ is preferred in the patterns 
Hla and HIb. The preferences of parameters we have got from the scan results can be justified by analysing 
the shapes of the s curve and the (/> curve and the intersections between the two curves in different patterns. 
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We also note that all the patterns prefer large quartic scalar couplings and moderate mixing angle between 
the Higgs and the scalar. 

Finally we combine the constraints from strong first order phase transition and the experimental limits 
on the S, T parameters, Higgs coupling measurements, and direct LHC searches. We found that there is still 
a signihcant amount of parameter region for the fermion mass and couplings to satisfy all the constraints, 
and have adequate CP violation strength to realize the baryon asymmetry at the same time. The new scalar 
with mass around 500 - 1100 GeV and mixing angle sint/? around 0.25 — 0.42 are still allowed and favored by 
the strong first order phase transition. The future Higgs coupling measurements and the Higgs boson pair 
production cross section will be able to further explore the allowed parameter space. 
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A Details of The Effective Potential 

A.l The Effective Potential in the On-shell Scheme 

In this appendix, we first review how the effective potential is written in on-shell scheme. Then we extend 
the discussion on SM Higgs boson in Ref. [45,46] to the Higgs-new scalar mixing case. 

We start from the one-loop effective potential in the Landau gauge, with the dimensional-regularization 
applied. 


Veff(T=0) = Ftree + Hcu^ 
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where the super-trace is taken among all the dynamical fields </? that have s) dependent masses M^. The 
UV divergent piece are defined as 
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The UV divergence has to be absorbed by the counterterms. We introduce the following counterterms 

AV = A{(j)^ — v^) + B{cj)^ — -I- C{s — u) -I- D{s — u)^ + E{(j)^ — v‘^){s — u) 

-\-F{a — uf' + G{s — u)'^ + — v‘^){s — u)^. (A.3) 


The renormalization conditions are needed to hx the above counterterms. In the MS renormalization scheme, 

the renormalization conditions consists in subtracting the term proportional to — 7e + logdyr in the 

^ 2 
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regularized potential. We will choose the on-shell renormalization scheme. The effective potential can be 
expanded using the one-particle irreducible (IPI) Green function at zero external momentum: 


yeff($) = -^($-$vAc)”r(”)(p = 0), 

n—0 


(A.4) 


where $ denotes the scalar fields (</>,«) and $vac denotes the vacuum {v,u). Therefore, we define the 
renormalized mass of the scalar field as the negative inverse propagator at zero momentum 


m 2 = -r(2)(p = 0) 
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Of course, we could also define the renormalized couplings as the four-point IPI Green function P^^^. However, 
the renormalization conditions on the couplings are not unique. Since we want to keep our discussion as 
general as possible, we only impose the tadpole conditions and mass conditions, as follows 
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(A.6) 


The five conditions fixes the tree level VEVs and scalar masses to be the physical ones, regardless of the 
couplings. As all the renormalization conditions are evaluated at point (j) = v and s = u, the only relevant 
variables here are only A, B,C, D, E, which can be uniquely determined, while the other three are totally 
arbitrary 

Although we can solve for these 5 coefficients using the 5 equations, but we found an easy way to do it 
systematically. The trick is to make use of the following function [45]: 
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where the mass matrices are supposed to be diagonal, in the basis of mass eigenstates, and M^ 

are evaluated at the vacuum point 4> = v, s = u. It’s easy to verify that the function satisfies all the 

5 renormalization conditions that we apply, regardless of the details of the mass matrices®. 

One may use the following counterterms 


Ay = As + Bs^ + Cs^ + Ds^ + Ecf)^ + F(f>^s + 


(A.7) 


However, in this parametrization, if we only apply the five conditions, it is not enough to determine certain counterterms. 
One has to use three renormalization conditions on the couplings to determine them uniquely. Due to the arbitrariness on the 

renormalization conditions on the couplings, the counterterms could be quite different. 

^ By adding higher powers of with appropriate coefficients, one can construct functions that satisfy higher- 

order on-shell conditions. For instance, if we want to have all the Lagrangian couplings to be the on-shell values, we simply 
add terms up to the 4th power of to it: 
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The problem is that whether it can be achieved from the original Coleman-Weinberg potential through 
adding counterterms like Eq. A.3. The answer is, luckily, yes for theories without mixing particles, but is no 
for the model we are dealing with, where both scalar sector and top quark sector may have large mixing. 
To see this, let’s take the difference 


AE = VSl - Vcw = 
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For mixed fields, though the mass matrix elements in gauge eigenbasis are usually polynomials of the scalar 
fields, which are allowed in the counterterms, in the mass eigenbasis they are typically irrational expressions. 
Specifically, TrM,^ ^ and TrM^ ^ can still be expressed in terms of the coefficients of the characteristic poly¬ 
nomial of the mass matrix, and hence are polynomial of the scalar fields, but TrDM^ ^ and TrDM^ ^ are 
not, where D is a diagonal matrix with non-degenerate eigenvalues. The two terms in the above expressions 
are exactly in this form. 

What we do is to expand AE at the vacuum point like in Eq. A.3, and truncate the expression at 
the order as we like. For instance, we can retain the terms of A' , B' , E', F', H' , and throw away all the 
other terms, so that the 5 renormalization conditions are still satisfied. The coefficients we get this way 
are unambiguous, which must be the same as what people get by any other methods. In our calculation, 
we retained all the 8 terms that are allowed in the counterterms, thus recover the full form of AE, and 
throw away the higher order terms. The coefficients C',D',G', however, are ambiguous, which depend on 
the additional renormalization conditions that people may add to the scheme. 

Suppose that after truncation, we get AE out of AE, therefore the final on-shell potential is 


in which the coefficients are 
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A.2 Goldstone Infrared Divergence 

There are two problems for the above potential even before the truncation. First, by dehnition, the potential 
is defined at scale = 0, and the second derivatives don’t give pole masses but renormalized masses at scale 
/i = 0. Second, in the Goldstone contribution to the Coleman-Weinberg potential, there is IR divergence 
from liimQ{v, u), since the field-dependent mass of the Goldstone boson 

mc(0, s) = - v^) + - u^) + p.s 4 >{s - m), (A.13) 
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is zero at the vacuum point. In this section, we will show that these two effects cancel each other, according 
to the discussion in [19]. 

Let’s find out the relation between the zero momentum parameters appearing in the effective potential 
and the physical observables that we need in an OS scheme. In general, we have the vertex functions 

+ TL{pi) + (A.14) 


where is the tree-level renormalized coupling, or the inverse propagator in the case of n = 2. is the 
loop contributions, which depends on the external momenta. comes from the counter terms defined in 
scheme X. 

As we are working in the OS scheme, we have 

,4>n\Pi) ='^phy+'^L{Pi)+'^°f (A.15) 

while 

rL(O5) + ro/ = 0 (A. 16) 

where the OS inside the parenthesis indicates the on-shell momenta, instead of on-shell scheme, as is 
scheme independent. Thus we have the = 0 values for the vertex functions: 

r(01, </)2, . . . , = 0) = Tj,hy + TLipf = 0) - rL(05) = Vphy - ATl (A.17) 


where AT^ is defined to be the difference between on-shell loop contribution and zero-momenta loop con¬ 
tributions. 

In this spirit, the renormalization conditions Eq. A.6 should be modified as 
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where S is the loop contribution to the mass matrix. Tadpole term is not changed because the tadpole loop 
does not depend on external momentum. For couplings, we have conditions like 
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p—V,S — U 

Note that, unlike the case for masses, the Coleman-Weinberg potential does have contributions to the 
couplings (unless we use the more complete form Eq. A. 9). 

One may notice that the off-diagonal element of E is not well-defined in OS scheme. In addition, some 
couplings like also don’t have natural on-shell definition. Here we assume that the AT^’s are not sensitive 
to the tricky details of how we define the OS renormalization conditions. Here we only focus on the IR 
divergence from the Goldstone loops, for which we use = ra^{v,u) as an IR regulator, and choose a 
convenient but inexact form of the IR-finite part. 

Thus we only retain the IR divergent terms, and replace all the parameters for physical quantities we 
need in an OS scheme. These terms are 
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while the Inm^ order IR divergences are: 
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where lnm|^ is for the IR finite contributions which, as already explained, are not exact, but errors are 
negligible. The power-law IR divergences in the couplings will be cancelled by the Coleman-Weinberg 
contributions (such as the last two term in Eq. A.9, which we don’t prove explicitly here. 

Plugging them into Eq. A.20, together with the knowledge of Eq. A. 13, we get 
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which neatly cancels the IR divergence of the Coleman-Weinberg term. Thus we only need to replace the 
Goldstone pole mass in the logarithm in Eq. A .8 with the higgs mass to fix the IR divergence problem. 
Again, there are corrections and small subtleties from momentum-shift effect, such as the AE contribution 
to the counterterms which are computable but we neglected. 


A.3 Field-dependent Masses 

In effective potential, the particles running in the loop are the particles in the model with the following 
degrees of freedom in the Landau gauge: 

nw = 6, nz = 3, = 3, Uh = ns = 1, fft = —12, nr = —12. (A.28) 

The field-dependent masses of the top quark, gauge bosons and Goldstone bosons at zero temperature are 
given by 

^ X, (A.29) 

+ + fis4,s. (A.30) 

The field-dependent masses of the scalars h and S are obtained as 

\ s) + rnl,{(l>, s)) T s) - s)) -h 4m4_^(0, s), (A.3I) 
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where the field-dependent quantities are 
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The field-dependent masses of the top quark and heavy vector-like top quark T are obtained as 

(A.35) 

where the field-dependent quantities are 
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(A.38) 


The finite-temperature potential needs to be corrected by the thermal field-dependent masses. The ther¬ 
mal field-dependent masses is calculated by adding the Debye masses, calculated from the the quadratically 
divergent bubbles and Daisy resummation. This leads to a shift of the bosonic field-dependent masses 
mf {(j), s) to the thermal field-dependent masses (Debye masses) 

w-(</>, s, T) = s) -b n,(^, s, T), (A.39) 


where s,T) is the self-energy of the bosonic field i in the IR limit. In particular, the longitudinal and 

transversal polarizations of the gauge bosons have to be taken into account separately: only the longitudinal 
components get a thermal mass correction and the transversal ones will not. Since the ring diagrams will 
only contribute significantly at high-temperature, only the zero-mode of the Matsubara frequency behave as 
a massless degree of freedom and generate IR-divergences at high-temperature, while other modes lead to 
subdominant contributions. For the SM bosonic contributions, the gauge boson thermal self energy is 
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For the new scalar bosons, 



M^{cj,,s,T) = M^icl,,s)+h 

\0 CsJ 

(A.44) 

where 

^0 , , 36f2 + ^'2 yf 

= 2 + 24 + 16 + 4 ’ 

(A.45) 
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(A.46) 

After diagonalization. 

we obtain the Deybe squared mass s,T). 

For the SM bosonic particles, we 

obtain 

^Wl = IIwt = 0, 

(A.47) 


Hz, = ^{9^+9'^)T^ Bz, =0, 
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(A.48) 


n,r = C^T^. 
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